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1 .  INTRODUCTION 

Recent  progress  in  the  theory  of  spectroscopy  and  dynamics  in  laser- 
irradiated  adspecies-surface  systems  is  discussed.  This  article  is 
divided  into  three  main  sections.  Beginning  with  Section  2  following  this 
Introduction,  a  first-principles  approach  for  describing  vibrational 
excitation  and  relaxation  of  a  laser-excited  adatom  is  presented.  A 
tractable  approximation,  known  as  the  isomnesic  (constant -memory) 
approximation  is  introduced  in  order  to  numerically  deal  with  the  long¬ 
time  behavior  of  the  dynamics.  In  Section  3  surface-dressed  optical  Bloch 
equations  are  solved  for  the  resonance  fluorescence  spectrum  of  a  laser - 
driven  two- level  atom  near  a  flat  metal  surface  and  also  near  a  rough 
metal  surface  modeled  by  a  hemispheroid  on  a  perfectly  conducting  plane. 

In  Section  4  the  problem  of  laser-induced  period  chemical  vapor  deposition 
is  addressed,  and  a  theoretical  analysis  of  selected  aspects  of  an 
experiment  on  the  deposition  of  cadmium  onto  a  silica  surface  is  carried 
out,  where  the  cadmium  is  produced  by  laser- induced  dissociation  of 
dimethyl  cadmium  gas  above  the  surface. 

2.  VIBRATIONAL  EXCITATION  AND  RELAXATION 

We  describe  here  a  first-principles  approach,  which  is  divided  into 
four  parts.  Part  2.1  is  a  brief  introduction  to  the  Zwanzig  formalism  for 
dealing  with  non-Markovian  phenomena  and  approximations  which  make  the 
method  capable  of  describing  long-range  chemical  processes  with  memory 
effects.  Part  2.2  is  an  application  of  the  theory  to  a  fictitious  system 
chosen  to  best  illustrate  effects  peculiar  to  resonantly  laser-excited 
vibrational  manifolds,  especially  effects  resulting  from  selective  laser 
photochemistry.  In  Part  2.3  the  theory  is  applied  to  a  "real"  system  -- 
hydrated  atomic  hydrogen  on  the  KCl(OOl)  surface  --  and  the  results  are 
extended  to  chemical  reaction  times.  Finally,  in  Part  2.4  a  desorption 


channel  is  introduced,  and  important  trends  in  the  resulting  desorption 
rate  are  discussed. 


2.1.  BASIC  THEORY 

The  theoretical  basis  for  the  first-principles  approach  is  the 
Heisenberg  equation  of  motion 

p(t)  »  ( ill)  *  -iLp(t)  , 


where  H  is  the  Hamiltonian  operator,  L  is  the  corresponding  Liouvillian 
operator,  and  p(t)  is  the  density  operator  of  the  system.  The  Hamiltonian 
is  partitioned  into  zeroth-order  Hamiltonians  for  the  solid  and  the  adbond 
(the  adspecies  in  the  phonon-averaged  field  of  all  the  lattice  atoms)  plus 
a  perturbation  representing  the  adbond- lattice  interactions  and  another 
representing  the  adbond-laser  interaction.  Within  the  Born  approximation 
and  random-phase  assumption,  the  equation  of  motion  fyr  an  adatom-surface 
system  leads  to  the  generalized  master  equation  (GME) 


ps(t) 


f  dt'  t  £  Kss’(t"t,)ps‘(t,) 


]>  Ks,s(t-t')Ps(t')]  ,  (2) 


which  is  an  integrodif ferential  equation  of  the  Volterra  type.  Pg  is  the 
occupation  probability  of  the  vibrational  state  |s>  of  the  adbond,  and 
the  memory  kernels  K  ,  incorporate  the  phonon  (p)  and  the  radiation  (r) 
terms:  K„„,  *  Kgj?,  Far.a  Debye  model  of^t^e  solid,  we  have 

obtained  the  phonon  portions  Kgg,  in  closed  form,  and  find  them  to  be 
fairly  complicated  with  a  predominantly  damped  oscillatory  behavior. 
Numerical  solution^of  a  Volterrra-type  equation  with  a  nonmonotonic  kernel 
is  ^ery  difficult.  In  our  case  we  were  able  to  solve  it  only  for  t  <  5 
ps.  As  a  general  rule,  we  have  found  that  exact  numerical  solutions  of 
the  problem  of  vibration  laser-adbond- surface  energy  transfer  are  possible 
only  for  times  of  the  order  of  a  few  Debye  periods  (inverse  of  the  Debye 
frequency,  .  This  is  hardly  sufficient  to  describe  chemical  processes, 
so  we  have  been  studying  possible  techniques  for  overcoming  this  "non- 
Markovian  bottleneck." 

One  such  technique  is  based  on  a  constant -memory  approximation,  i.e., 
isomnesic  approximation  (IA),  to  the  memory  kernel.  Briefly,  the  short 
temporal  range  of  lCg,(t)  suggests  a  delta-function  representation  if  the 
long-time  behavior  {accurate  beyond  this  temporal  range)  is  of  interest. 
Thus, 

Kss,(t)  *  Qss,6(t)  +  K^!  .  (3) 

The  second  term  on  the  right-hand  side  is  time- independent.  If  it  were  a 
delta  function  also,  we  would  obtain  the, Markovian  approximation  (MA). 
However,  because  of  the  constant  terms  Kgg,  we  obtain  instead  a 
nonconvolution-type  integrodif ferential  equation 


ps(t) 


\  Kss'  {  dt'  ps ' ^ +  \  nss,ps'^ 


which  reduces  to  the  following  second-order  differential  equation  (we  note 
that  the  Markovian  master  equation  is  first  order): 

hMml  Kss!  ps'(t)  +  Inss'fs'(t)  ■  (5) 

S'  S' 

We  have  found  a  convenient  method  for  solving  t^is  based  on  Laplace 
transforms  and  the  Heaviside  expansion  theorem. 

2.2.  SELECTIVE  PHOTOEXCITATION  OP  A  PROTOTYPE  SYSTEM 

As  a  prototype  system,  we  choose  a  mass  16  species  on  Ge  which 
develops  six  bound  states  (S  =  0,1, 2, 3, 4, 5)  in  the  adbond  potential.  We 
also  choose  a  situation  where  |u)jq|  =  uj^  and  Jw.^l  *  where  «>£<.,  is  the 
transition  frequency  between  levels  S  and  S'  ana  is  the  laser 
frequency,  so  that  the  laser  does  not  cause  direct  transitions  between 
phonon  states.  The  resulting  population  profiles  are  shown  in  Figure  1 
for  S  *  0,  1  and  3  along  with  results  of  the  exact  numerical  solution  of 
the  GME,  Eq.  (2).  As  expected,  there  is  a  loss  in  fine  detail  with  the 
1A,  but  the  average  behavior  is  surprisingly  well  reproduced.  Thus,  for 
example,  the  crossing  point  where  P^(t)  =  P^(t)  for  the  first  time  is  at  t 

*  1.7  ps  in  the  IA  as  compared  to  t  *  2.2  ps  in  the  exact  numerical 

solution.  The  relative  behavior  of  P^,  P^  and  P^  is  also  close  to  that 
obtained  numerically. 

Substantial  chemical  information  is  contained  in  the  average  adbond 
energy 

E(t)  -  £  Ps(t)Es  ,  (6) 

S 

where 

HS|S>  *  ES|S>  .  (7) 

Figure  2  depicts  E(t)  obtained  by  solving  the  exact  and  isomnesic  GME.  We 
notice  that  the  energy  does  not  change  monotonically,  but  consists  of 
alternating  influx  and  efflux.  Figures  1  and  2  together  display  a 
feedback  mechanism.  Depletion  of  level  S  3  0  corresponds  to  energy 
transfer  out  of  the  phonon  field,  and  we  see  such  a  depletion  for  t  <  5 
ps.  This  feedback  is  an  ^mportant  microscopic  effect  suggested  in  our 
earlier  model  calculation  but  often  missed  in  average  statistical 
theories  which  examine  only  the  long-time  steady-state  behavior  of  many- 
body  systems.  Figures  3  and  A,  which  show  P<,(t)  and  E(t)  using  the  IA  for 
20  ps,  further  clarify  the  situation.  The  average  energy  content  of  the 
adbond  is  periodic,  with  an  initial  fast  rise  followed  by  a  fall  and 
subsequent  smaller-amplitude  damped  oscillations.  Between  0  and  A. 8  ps, 
Pg(t)  decreases  while  E(t)  increases;  energy  is  transferred  from  the  solid 
(and  the  laser,  of  course)  to  the  adbond.  However,  subsequently  between 
A. 8  and  7.7  ps,  PQ( t )  increases  while  E(t)  decreases,  i.e.,  the  phonon 
field  is  replenished  at  the  expense  of  the  adbond  energy.  But  then  again, 
the  laser  causes  sufficient  depletion  of  P.  for  phonon  energy  to  get 
transferred  back  to  the  adbond,  causing  a  depletion  of  PQ,  which  is  energy 
feedback.  The  solutions  of  a  Markovian  version  of  this  formulation 


E(t)  (m 


EXACT  vs.  ISOMNESIC 


FIGURE  1.  Probability  profiles  Pg(t)  for  levels  S  *  0,  1  and  3  of  the 
adbond.  Solid  lines:  exact  numerical  solutions  of  the  GME,  Eq.  (2); 
broken  lines:  closed-form  solutions  with  the  isomnesic  approximation. 


FIGURE  2.  Average  adbond  energy  E(t)  obtained  via  the  exact  numerical 
solution  of  Eq.  (1)  (solid  line)  and  closed-form  solution  with  the  IA 
(broken  line). 


(obtained  by  replacing  Kgg,  in  Eq.  (3)  by  KSg,x6(t),  where  x  is  a  time- 
range  parameter)  display  no  such  feedback  mechanism  (broken  lines  in 
Figures  3  and  4). 

The  isomnesic  approximation  represents  a  significant  improvement  over 
the  Markovian  approximation  and  has  definite  advantages  over  the  exact 
numerical  treatment,  since  the  IA  is  an  analytic  (nonnumerical)  theory  and 
is  not  subject  to  the  time  constraints  inherent  in  the  exact  first- 
principles  treatment.  We  are  thus  able  to  follow  the  time  evolution  of 
the  adbond  probability  distribution  for  arbitrarily  long  times  and  to 
study  the  steady  state.  The  Markovian  limit  is  studied  for  the  purpose  of 
establishing  a  link  between  our  new  theory  and  the  majority  of  past  work 
on  the  chemical  dynamics  of  many-body  systems  far  from  equilibrium.  We 
find  that  such  a  (Markovian)  treatment  is  inherently  arbitrary  and  leads 
to  qualitatively  different  results  depending  on  the  choice  of  the  time- 
range  parameter  t  for  the  laser  radiation.  Approach  to  equilibrium  in 
the  Markovian  approximation  is  monotonic  (a  superposition  of  exponential 
changes),  while  the  exact  and  IA  treatments  give  rise  to  oscillatory 
behavior  for  the  population  distribution  of  the  adbond.  However,  we  do 
find  that  the  steady  state  is  the  same  for  the  IA  and  the  MA.  The 
important  point  here  is  that  the  IA  is  essentially  as  convenient  and 
simple  to  work  with  as  the  MA,  but  gives  results  which  are  very  faithful 
to  the  exact  first-principles  description  of  the  phenomenon.  We  now 
proceed  to  applications  of  the  above  approach  involving  the  IA. 

2.3.  LONG-TIME  EFFECTS  OF  SELECTIVE  PHOTOEXCITATION 


Recent  studies  have  suggested  that  extremely  high  laser  powers  would 
be  necessary  to  effectively  desorb  adspecies  from  a  solid^surface  via  a 
selective  absorption  mechanism.  However,  our  recent  workJ  has  shown  that 
a  very  small  degree  of  detuning  between  the  laser  and  energy  levels  of  the 
system  excludes  the  possibility  of  long-term  energy  absorption.  The 
reason  for  this  behavior  is  seen  to  be  the  oscillatory  nature  of  the  laser 
energy  absorption  "rate"  (actually,  the  memory  kernel)  for  the  adbond, 
with  a  frequency  related  to  the  detunjgg,  ty^jcaljy  a  small  percentage  of 
the  laser  frequency  itself,  namely  10A  -  l(j  s  .  Substantial 

cancellation  in  the  overall  absorption  results  and  a  very  high  laser  power 
would  be  needed  for  the  absorption  to  compete  with  phonon  relaxation 
effects.  In  order  to  emphasize  the  special  role  of  the  coherence  and 
monochromoticity  of  lasers  in  selective  photoexcitation,  a  system  with 
essentially  no  detuning  has  been  chosen  for  the  application  of  the 
isomnesic  approximation  to  long  time  IR  laser  pumping.  q 

Our  formalism  has  been  applied  to  the  system  H(H20)/KC1(001),  where 
the  primary  reasons  for  choosing  this  are  the  shallow  adsorptio^potential 
well  and  the  small  number  of  bound  states  (seven)  it  generates.  Figure 
5  is  a  composite  display  of  the  function  E(t)  for  six  ranges  o£  time 
differing  by  factors  of  10  and  for  a  laser  intensity  of  1  W/cm  .  It  is 
seen  that  in  going  from  a  picosecond  time  range  to  a  microsecond  time 
range,  the  behavior  of  the  system  changes  dramatically.  The  early 
behavior  is  dominated  by  fast  energy  transfer  to  phonon  modes,  completely 
overwhelming  the  slow  laser  pumping.  The  monotonic  decay  is 
representative  of  relaxation  phenomena.  For  the  100  ps  range,  a  false 
steady  state  seems  to  have  been  reached,  but  is  in  fact  only  an  artifact 
of  the  large  difference  between  the  amplitudes  of  the  laser  and  phonon 
kernels.  The  adbond  energy  starts  to  increase  in  the  rts  range,  again 
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FIGURE  5.  E(t)  obtained  by  using  the  isomnesic  approximation  for  time 
scales  ranging  from  ps  to  ps  for  a  laser  intensity  of  1  W/cm  . 
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FIGURE  6.  Same  as  Figure  5  for  a  laser  intensity  of  25  W/cm 


suggesting  a  possible  steady  state.  However,  this  is  a  transitional 
period  during  which  the  relative  importance  of  the  phonon  terms  and  the 


laser  terms  starts  to  change.  The  letter,  in  fact,  originate  from  a 
probability  profile  of  the  form  (Ae  cosgt  -  |J  sinfit),  where  a  +  ip  are 


complex  roots  of  det  M(s)  =  0,  and  A  +  ip  are  complex  coefficients 
obtained  from  M(s)  and  ^(s),  where  the  matrix  M(s)  and  the  vector  |(s)  are 
related  to  the  Laplace  transform  of  the  isomnesic  master  equation. 

The  phonon  terms  come  from  real  roots  and  are  therefore 
nonoscillatory .  The  magnitudes  of  the  exponents  a  in  the  laser  part  of 
the  solution  are  much  smaller  than  corresponding  exponents  in  the  phonon 
term.  This  is  responsible  for  the  radically  different  nature  of  E(t)  in 
the  different  time  regimes.  Thus,  the  seemingly  monotonic  rise  in  E(t) 
apparent  in  the  10  ns  regime  is  seen  to  be  only  the  early  segment  of  a 


cycle  which  becomes  evident  in  the  100  ns  results  and  beyond.  The 

c „  r  _  l  ^  ^  v  i  „  u  — . 


frequency  of  this  oscillation  is  seen  to  be  Kq  ,  where  Kq  is  the 
amplitude  of  the  laser  kernel  Kggt,  and  is  essentially  the  Rabi  frequency 
for  the  pair  of  levels  resonant  with  the  laser. 


While  the  frequency  of  the  long-term  oscillations  of  P(t)  and  E(t) 
depends  only  on  the  laser  term,  being  proportional  to  the  square  root  of 
the  intensity,  the  rate  of  onset  of  this  behavior  depends  on  a  which 
varies  linearly  as  Kq  and  inversely  as  the  amplitude  of  the  ghonon 
kernel.  Thus,  increasing  the  laser  intensity  from  1  to  25  W/cin  advances 
the  onset  time  for  the  oscillatory  state  while  increasing  the  frequency  by 
a  factor  of  5  (as  displayed  in  Figure  6). 

Physically  the  results  of  Figures  5  and  6  can  be^nterpretec^in  terms 
of  desorption  rates.  Thus  the  time  required  for  | A | e  and  |p|e  to 

become  large  enough,  so  that  the  amplitude  of  E(t)  in  the  oscillatory 
regime  is  substantially  larger  than  the  (t  *  0)  equilibrium  Boltzmann 


value,  is  a  reasonable  measure  of  the  desorption  time  to  within  an  order 

_  C _ z  a, 1  _  Tr|_  f  _  .  • z  a.  1  i.  .  J  1 _  _  A  T3 _ _ 


of  magnitude.  This  time  is  estimated  by 


XD  =  a 


For  our  case, 


0.5  ^  5  (10*  I)s' 


where  w  is  the  laser  intensity  in  W/cm  .  Thus,  for  a  1  W/cm 


lasgr,  a 
10  s  on 


hydrogen  atom  will  stay  adsorbed  on  the  KCHH^O)  surface  for  -10  s  on 
the  average  before  desorbing,  or  at  least  becoming  very  highly  excited. 


The  data  mentioned  here  is  for  T  *  150  K. 

It  must  be  pointed  out  that  a  number  of  mechanisms  whereby  the 
desorption  cross  section  could  be  modified  are  not  included.  Spontaneous 
decay,  phase  relaxation,  phonon-phonon  interaction  and  electron-hole  pair 
creation  are  some  of  the  ones  being  considered  in  extensions  of  this  work. 
As  a  result,  the  sharp  transitions  assumed  here  will  be  broadened  by  a 
detuning-type  effect,  and  estimates  of  the  desorption  rate  presented  above 
will  have  to  be  modified.  However,  the  selectivity  represented  by  the 


dominant  T^  energy-transfer  process  will  be  retained. 


2.4.  DESORPTION 


Desorption  is  the  most  commonly  studied  phenomey^n^f  laser- 
stimulated  surface  processes,  and  recent  experiments1  have  shown  a 
strong  dependence  of  the  desorption  cross  section  on  the  laser  frequency. 
While  energy  transfer  during  IR-laser  stimulated  surface  processes  has 
been  discussed  for  the  above  work  in  the  absence  of  a  desorption  channel, 
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selective  absorption  has  been  demonstrated  for  vanishing  mismatch  between 
the  laser  and  adbond  frequencies.  The  time -dependent  energy  profiles  for 
the  adbonds  showed  an  initial  strong  increase  followed  by  a  decrease 
(transfer  to  bulk  phonons)  and,  subsequently,  another  rise  due  to  feedback 
from  the  phonons.  We  now  present  the  results  of  an  isomnesic  formulation 
including  desorption  in  the  form  of  a  tran^tion  from  an  excited  state  of 
the  adsorption  potential  to  the  continuum. 

Let  us  consider  a  three-level  system  --  |S>  =  |0>,  |l>  and  |3>  --  one 
of  whose  transitions  is  resonant  with  the  laser,  and  another  equals  the 
phonon  Debye  frequency.  In  the  absence  of  desorption,  application  of  the 
Laplace  transform  method  to  the  GME  leads  to  a  third-order  polynomial 
equation.  With  the  inclusion  of  a  depletive  desorption  channel,  the  GME 
reduces  to  a  quartic.  As  expected,  the  total  population  for  the  adbond  is 
depleted  monotonically ,  and  the  desorption  probability  Pn(t)  grows 
correspondingly  to  a  maximum  value  PT(0).  The  average  aabond  energy  E(t) 
increases  with  time,  but  because  of  the  competition  between  laser  pumping 
and  desorption,  the  detailed  behavior  of  the  energy  profile  is  not  as 
simple  as  that  of  P^(t).  Depending  on  the  magnitude  of  the  bound-to- 
continuum  transition  rate  R,  E(t)  has  a  monotonic  (nonoscillatory) , 
sigmoid  profile  which  is  is  an  oscillatory  profile  akin  to _^e 
non^esorptive  behavior,  or  a  transitional  profile.  For  10  ps  <  R  <  1 
ps  ,  both  the  cessation  of  slow  oscillatory  behavior  and  the  onset  of 
fast  sigmoid  behavior  take  place  at  higher  laser  powers  for  larger  values 
of  R.  The  trends  observed  are  important  from  the  point  of  view  of 
selective  laser  photochemistry.  Thus,  for  a  given  value  of  R,  laser 
powers  falling  in  the  oscillatory  region  will  cause  a  very  slow 
accumulation  of  vibrational  energy  in  the  adbond,  while  laser  powers 
falling  in  the  transitional  or  sigmoid  regions  will  excite  the  adbond 
quickly  and  cause  almost  instantaneous  desorption. 

3.  RESONANCE  FLUORESCENCE  OF  A  TWO-LEVEL  ATOM  NEAR  A  METAL  SURFACE 

The  lifetime  of  an  excited  molecule  has  experimentally  been  |gund  to 
vary  dramatically  as  a  function  of  distance  from  a  metal  surface.  A 
number  of  researchers  have  examined  this  effect  from  the  viewpoint  of 
ref lected-f ield  theory,  where  the  basic  calculation  is  concerned  with  the 
interaction  between  an  excited  molecule  and  its  own  reflected  radiation 
field.  This^h^gry  generally  provides  good  agreement  with 
experiments. 

When  an  adatom  is  driven  by  strong  resonant,  driving  coherent  field, 
it  creates  for  the  atom  or  molecule  an  environment  where  the  probability 
of  stimulated  emission  can  exceed  that  of  spontaneous  emission.  Under 
this  condition,  dynamic  ac-Stark  splitting  and  nutational  oscillation  of 
the  emitted  light  intensity  become  important  parts  of  the  laser-driven 
processes,  such  that  interesting  "resorj^nce  fluorescence"  and  other 
nonlinear  optical  phenomena  can  occur.  recently  derived  a  set 

of  surface-dressed  optical  Bloch  equations,  by  which  we  can  examine 

the  dynamics  and  relaxation  of  an  adatom  near  a  metal  surface.  By  solving 
these  equations,  we  have  been  able  to  evaluate  the  resonance  fluorescence 
spectrum  of  a  two-level  atom  near  a  metal  surface,  taking  the  following 
factors  into  account: 

(i)  the  reflected  electromagnetic  field  effects  due  to  the 
presence  of  the  interface 
(ii)  collision  dephasing 
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(iii)  surface- induced  dephasing  due  to  reflected  photons 

(iv)  resonance  excitation  of  surface  plasmons 

(v)  random-phase  fluctuations  of  the  laser. 

Besides  these  studies,  we  have  been  extending  the  flat-surface  resonance 
fluorescence  calculations  to  the  rough  surface  case,  which  modeled  as  a 
hemispheroidal  profusion  on  a  perfectly-conducting  surface.  We  discuss 
our  progress  below. 

Let  us  begin  with  the  case  of  a  flat  metal  surface  and  the  set  of 
surface-dressed  optical  Bloch  equations  (SBE),  which  include  the  effect  of 
surface-reflected  photons,  i.^j .^ghotons  emitted  by  the  laser-driven  atom 
and  reflected  by  the  surface.  ’  Within  the  rotating-wave  approximation 
the  SBE  take  the  form 
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where  ft  =  |2><2|  -  1 1 X 1 1  is  the  population  inversion  of  t^e  atom,  S.. 
product  of  exp(iuct)  and  the  transition  operator  |i><j|,  fi-(t)  is  the'' 
time-dependent  Rafei  frequency,  A  =  tu-j  -  is  the  detuning  of  the  laser 
(co^)  with  respect  to  the  two- level  atom  (a^, ) »  and  the  dephasing  rate 
constant  y^  the  sum  of  the  radiative  depnasing  a°d  the  surface- 
induced  dephasing  y  .  The  driving  and  reflected  fields  are  treated 
semiclassically,  an§  we  assume  the  atom-surface  distance  to  be  large  (>,  30 
nm).  The  relaxation  time  included  in  the  Drude  model  for  the  complex 
dielectric  constant  of  the  metal  medium  represents  the  dissipation  of 
electron  gas  which,  together  with  surface  plasmon  resonances,  influences 
the  surface  reflectivity  and  hence  the  behavior  of  the  reflected  field. 
Here  we  may  neglect  nonradiative  transfer  of  energy  from  the  excited  atom 
to  the  metal.  Effects  of  the  laser  bandwidth  are  included  by  means  of  a 
phase  diffusion  model  for  the  driving  field.  In  the  weak-field  or  large- 
detuning  limit,  the  power  spectrum  of  scattered  light  has  two  peaks:  one 
corresponding  to  Rayleigh  scattering  at  the  laser  frequency  and  the 
other  to  fluorescence  at  the  atomic  transition  frequency  .  For  a 
sufficiently  strong  driving  laser  field,  the  spectrum  exhibits  three 
peaks:  a  central  one  at  (Rayleigh  component),  a  left  one  at  u>_  = 

2w-  -  "  $  (three-photon  component)  and  a  right  one  at  u>+  =  +  6 

(fluorescence  component),  where  6  is  the  ac-Stark  shift. 

Results  for  the  peak  heights  H_,  and  H+  are  shown  in  Figure  7  for 
the  case  of  a  silver  surface.  The  key  feature,  which  is  a  unique 
behavior  due  to  the  surface,  is  that  for  certain  atom-surface  distances  H 
is  larger  than  H  ,  whereas  in  the  pure  gas-phase  resonance  fluorescence 
spectrum  H  is  always  less  than  H+  (due  to  molecular  collisions  for 
positive  detuning).  Also,  the  population  inversion  of  the  adatom  and  the 
resonance  fluorescence  spectrum,  as  well  as  the  surface- induced  phase- 
decay  constant  of  the  adatom,  show  strong^o^illatory  behavior  as  a 
function  of  the  adatom-surface  distance.  ’ 

For  a  rough  surface,  let  us  considgg  ^surface  protusion  modeled  by  a 
prolate  hemispheroid  on  top  of  a  plane.  ’  We  note  that^his  model  has 
been  shown  to  be  identical  to  a  full  spheroid  in  a  vacuum,  so  that  our 
calculations  can  also  be  used  for  an  ellipsoidal  cluster.  The  two-level 


FIGURE  7.  Atom-surface  distance  dependence  of  the  heights  of  the  three 
peaks  in  the  resonance  fluorescence  spectrum  for  the  case  of  a  silver 
surface,  where  H_,  Hq  and  H+  correspond  to  the  peak  heights  (in  arbitrary 
units)  of  the  three-photon,  incoherent  Rayleigh  and  fluorescence 
components,  respectively.  D  ■  2kd  is  the  reduced  distance,  where  k  is  the 
wavenumber  and  d  is  the  atom-surface  distance.  The  Drude  model  is  used  for 
the  dielectric  function  of  the  metal,  while  the  dielectric  function  of  the 
gas  medium  is  set  equal  to  one.  The  detuning  -  uu^  is  1  A  and  the  Rabi 
frequency  is  10  A,  in  the  unit  of  Einstein's  A  coefficient.  The  solid 
curves  are  for  the  induced  transition  dipole  of  the  atom  oriented 
perpendicular  to  the  surface,  and  the  dashed  curves  are  for  the  parallel 


adatom,  which  is  located  at  a  distance  d  from  the  top  of  the  hemispheroid, 
is  driven  by  a  laser.  The  prolated  spheroidal  coordinates  (£»n>$)  are 
used  to  calcj^late^t^e^ref lected  field.  We  define  ^  *  (a  +  d)/f,  *  a/f 

and  f  =  (a  -  bz)  ,  where  a  and  b  are  the  semi-major  and  semi-minor 
axes  of  the  hemispheroid.  The  reflected  field  at  the  position  of  the 
ad^gom  (transition  dipole)  in  the  near-field  approximation  can  be  written 


Er  -  *  ll  W<1>  + 
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where  Q  denotes  the  Legendre  function  of  the  second  kind  and  p  is  the 
inducedndipole  moment.  The  expansion  C  depends  on  the  laser 

amplitude  E^  and  parame^r^CQ  and  Based  on  ?he  surface-dressed 

optical  Bloch  equations  ’  appropriate  for  the  excitation  and 
dissipation  of  a  two-level  adatom  near  a  hemispheroid,  we  fi^g  the 
surface- induced  phase  relaxation  constant  to  be  given  by 


Ys  =  ( 2/yi )  Im(  F ) 


where 

(1-e)?  Q’(C.) 

F  =  W  IeQ1(5(J)  -  C1Q|(^0)+  r  '  (12) 

Here  e  is  the  complex  dielectric  constant,  and  T  is  a  somewhat  complicated 
function  of  £  and  e. 

A  sharp  resonance  enhancement  of  the  adatom-hemispheroid  interaction, 
through  the  reflected  field  at  the  atomic  site,  occurs  when  the  specific 
shape  of  the  prolate  hemispheroid  corresponds  to  a  resonant  excitation  of 
plasmons  (see  Figure  8).  The  time  oscillation  of  the  level  population 
decreases  as  the  shape  of  the  hemispheroid  approaches  the  plasmon 
resonance.  We  also  find  that  the  strong-field  three-peak  fluorescence 
spectrum  is  strongly  influenced  by  the  roughness  of  the  surface.  The 
resonance  excitation  of  the  plasmon  in  the  hemispheroid  remarkably 
enhances  the  adatom-surface  coupling,  such  that  the  dephasing  processes 
broaden  the  linewidths  of  the  spectrum.  In  the  small  detuning  case,  the 
spectrum  has  a  distinctive  three-peak  nature,  where  the  three-photon  side 
peak  has  a  measurable  height  which  is  almost  comparable  to  the 
fluorescence  side  peak.  The  plasmon  and  ref lected-f ield  broadening 
influence  equally  the  three  spectral  peaks.  In  the  large  detuning  case, 
the  height  of  the  three-photon  peak  is  decreased,  and  the  three-peak 
spectrum  is  transformed  to  the  weak-field  two-peak  structure. 

A.  USER- INDUCED  PERIODIC  CHEMICAL  VAPOR  DEPOSITION 

We  have  been  engaged  in  theoretical  research  to  explai^ laser- induced 
periodic  vapor  deposition  recently  observed  experimentally.  While  some 
pieces  of  the  puzzle  remain  unsolved,  we  think  that  we  now  have  a 
sufficient  underspending  of  the  phenomena  to  warrant  a  more  general 
discussion. 

Briefly,  the  experiment  which  we  are  attempting  to  model  is  as 
follows.  A  dilute  organometallic  gas  (CdtClL,^  in  this  case)  is  placed 
above  a  silicon  or  silicon  oxide  surface.  The  system  is  illuminated  with 
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FIGURE  8.  Rough  surface- induced  phase-decay  constant  y  as  a  function  of 
the  semi-minor  axis  b  of  the  surface  protusion,  where  tfie  semi-major  axis 
is  fixed  at  100  A.  y  is  in  the  unit  of  Einstein's  spontaneous  decay 
constant  A. 
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FIGURE  9.  Graph  of  the  plasmon  field  intensity  as  a  function  of  grating 
wavenumber.  The  grating  height  is  25  nm,  and  15  terms  are  kept  when 
evaluating  the  modified  Bessel  functions. 
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a  UV  laser  (257  nm)  at  low  intensity  (-  10  W/cm  ).  This  dissociates  the 
organometallic  to  a  metal  radical  (in  this  case,  cadmium),  which  is 
deposited  on  the  surface.  Deposition  occurs  evenly  to  a  depth  of  about 
200  nm,  after  which  the  deposition  is  periodic,  with  ripples  developing  in 
the  p-polarized  direction  and  growing  to  an  added  height  of  about  100  nm. 
In  addition,  much  smaller  ripples  are  formed  in  the  s-polarized  direction. 
This  structure  is  very  distinct  and  readily  visible  under  a  scanning 
electron  micrograph.  ^ 

The  original  experimenters,  Brueck  and  Ehrlich, J  postulated  that  the 
laser  induces  a  plasmon  field  on  the  surface,  where  the  varying  field 
strength  leads  to  the  observed  periodic  deposition.  They  assumed  the 
deposition  rate  to  ^  proportional  to  the  plasmon  intensity,  and  using  a 
first-order  theory,  they  showed  that  the  plasmon  field  intensity  is 
proportional  to  the  grating  height.  It  is  therefore  easy  to  conclude  that 
the  grating  growth  is  exponential  for  shallow  gratings. 

We  accept  Brueck  and  Ehrlich's  hypothesis  as  the  starting  point  for 
our  own  work.  Their  discussion  raises  the  following  questions,  some  of 
which  we  shall  be  able  to  answer  in  this  article: 

1)  What  is  the  mechanism  by  which  the  plasmon  enhances  the 
deposition  process?  This  is  the  major  piece  of  the  puzzle  which  remains 
unsolved,  but  we  have  more  to  say  about  this  later. 

2)  What  is  the  sequence  of  events  in  this  experiment?  Is  the 
organometallic  molecule  dissociated  in  the  gas  phase,  with  deposition  of 
the  cadmium  atom  occuring  later,  or  is  the  molecule  physisorbed  onto  the 
surface,  with  dissociation  occuring  afterwards,  perhaps  as  part  of  the 
chemisorption  reaction  of  cadmium  with  the  surface? 

3)  What  is  the  dynamics  of  grating  growth,  and  why  does  the  grating 
stop  growing?  This  last  question  is  the  one  about  which  we  now  have  the 
clearest  understanding.  We  continue  to  work  on  the  first  two  problems  and 
shall  discuss  them  in  due  course. 

By  accepting  Brueck  and  Ehrlich's  hypothesis  that  the  grating  growth 
is  proportional  to  the  strength  of  the  plasmon  electric  field,  then  in 
order  to  derive  a  more  precise  expression  for  the  dynamics,  we  must  solve 
for  the  electric  field  above  a  rough  surface.  It  is  possible  to  do  this 
numerically  as  accurately  as  required,  but  in  order  to  gain  some  physical 
insight  into  the  problem,  we  make  two  key  approximation  which  permit  us 
to  write  an  analytical  expression  for  the  growth  rate: 

1)  We  assume  the  Rayleigh  hypothesis,  namely  that  the  solutions  to 
the  homogenous  Helmholtz  equation,  which  are  exact  above  the  selvedge 
region,  are  also  valid  within  the  selvedge  region.  This  holds  strictly 
only  in  the  limit  as  the  grating  is  infinitely  shallow.  How  valid  this 
approximation  is  depends  on  the  precise  nature  of  the  mechanism  by  which 
the  deposition  occurs.  We  note  that  the  gas  pressure  is  about  1  torr, 
which  means  that  the  mean  free  path  of  the  molecules  is  on  the  order  of  a 
millimeter,  much  longer  than  any  grating  height  measured  in  the 
experiment.  Thus  if  the  relevant  forces  involved  are  long-range,  most  of 
the  "action"  takes  place  well  above  the  selvedge  region  where  our  solution 
is  exact,  and  we  can  use  the  Rayleigh  hypothesis  with  impunity.  On  the 
other  hand,  if  the  mechanism  is  short  range,  then  this  introduces  some 
uncertainty  into  our  result,  but  we  still  expect  to  be  able  to 
qualitatively  account  for  the  phenomena. 

2)  Our  second  and  more  important  approximation  is  the  m^jiner  in 
which  we  truncate  the  Fresnel  matrix.  Toigo,  Marvin  and  Celli  derived 
uncoupled  equations  for  the  electric  field  strength  above  a  grating  (given 


the  Rayleigh^yggthesis) ,  and  their  result  has  been  used  by  many  other 
researchers.  '  In  principle,  it  is  necessary  to  solve  an  infinite 
number  of  equations  for  an  infinite  number  of  unknowns,  the  unknowns  being 
the  amplifies  of  the  various  order  plasmons  and/or  Bragg  reflections. 
Maradudin  has  used  this  expression  to  derive  the  dispersion  relation  for 
a  shallow  grating.  He  did  this  by  1)  assuming  that  the  grating  is 
shallow,  and  therefore  only  terms  to  first  order  in  grating  depth  need  be 
retained,  and  2)  that  because  the  resonance  condition  is  met,  all  field 
strengths  are  small  relative  to  the  resonance  plasmon.  This  means  that  we 
can  neglect  all  other  reflections,  and  our  infinite  set  of  coupled 
equations  reduces  to  two  equations,  easily  soluble  analytically. 

Our  contribution  to  this  result  is  to  notice  that  Maradudin* s  first 
assumption  is  unnecessary.  If  the  resonance  condition  is  met,  then  the 
resonance  plasmon  will  be  orders  of  magnitude  larger  than  the  other 
reflections,  and  hence  we  are  justified  in  keeping  all  orders  of  the 
grating  depth.  This  is  true,  of  course,  only  as  long  as  the  resonance 
condition  is  met,  but  it  implies  that  we  can  extend  Maradudin* s  work  to 
the  deep  grating  case. 

The  implications  of  our  approximation  are  best  illustrated  by 
studying  Figure  9,  where  the  x-axis  is  the  grating  wavenumber  and  the  y- 
axis  is  the  intensity  of  the  plasmon  field.  Throughout  we  are  insisting 
that  the  grating  depth  be  25  nm.  Thus  as  the  wavenumber  becomes  larger, 
the  grating  period  becomes  shorter,  and  the  relative  grating  height 
becomes  larger.  Hence  we  see  an  increase  in  the  plasmon  field  intensity, 
consistent  with  the  previous  result  that  the  plasmon  field  increases 
linearly  with  grating  height.  The  difficulty  here  is  that  this  increase 
appears  far  from  any  resonance.  Hence  our  approximation  means  that  the 
two  terms  that  we  are  keeping  are  on  the  same  order  of  magnitude  as  the 
terms  that  we  are  dropping.  This  leads  to  a  monotonically  increasing 
function  for  the  plasmon  field  strength,  when  in  reality  the  curve  would 
have  more  structure  and  would  tend  to  approach  some  asymptotic  value.  In 
summary,  for  grating  frequencies  off  resonance,  our  approximation  over¬ 
estimates  the  plasmon  field.  Later  we  shall  provide  some  measure  for  the 
degree  of  this  error. 

At  resonance,  on  the  other  hand,  the  situation  is  different.  Here 
the  two  terms  that  we  are  keeping  are  much  larger  than  any  others.  Thus 
our  approximation,  while  still  tending  to  over-estimate  the  field 
strength,  is  much  more  accurate.  Figure  10  illustrates  the  results  of  our 
method  for  various  grating  heights.  Note  the  structure  in  the  curves. for 
gratings  15  nm  and  20  nm  deep.  We  suggest  that  this  is  physically 
reasonable  since  the  resonance  is  very  large,  and  the  two  terms  we  keep 
are  quite  complicated  functions.  On  the  other  hand,  our  result  for  the  40 
nm  deep  grating  is  suspect  since  there  is  no  resonance  whatsoever.  In  all 
cases  we  have  kept  15  terms  in  evaluating  the  Bessel  functions  that  appear 
in  the  equations. 

Given  the  assumption  that  the  grating  growth  rate  is  proportional  to 
the  plasmon  field  intensity,  we  are  now  in  position  to  calculate  the 
growth  dynamics.  Two  problems  are  now  presented.  First,  until  we  have  a 
better  knowledge  of  the  mechanism  by  which  the  plasmon  increases  the 
growth,  it  is  impossible  to  calculate  the  proportionality  coefficient. 
Therefore,  in  our  calculation  shown  in  Figure  11,  we  have  simply  chosen 
the  coefficient  to  match  the  time-scale  of  Brueck  and  Ehrlich's 
experiment.  We  are  considering  several  possible  models  to  explain  this 
mechanism,  and  these  will  be  discussed  shortly. 


The  second  problem  is  to  determine  when  the  grating  stops  growing. 
This  happens  when  the  resonance  disappears  altogether,  and  gratings  of  all 
frequencies  grow  at  the  same  rate.  Thus  we  have  to  subtract  the 
"baseline"  from  the  plasmon  field  intensity,  namely  the  average  intensity 
of  all  other  plasmon  frequencies.  We  have  done  this  by  treating  the 
wavenumber  of  the  light  as  the  "base"  wavenumber.  This  value  is  the  left¬ 
most  point  of  the  curves  shown  in  Figure  10.  The  reason  for  choosing  this 
wavenumber  is  that  it  can  never  be  at  resonance  since  the  laser- induced 
plasmons  must  oscillate  at  a  wavenumber  greater  than  that  of  the  light. 
However,  our  approximation  over-estimates  the  intensity  of  these  off- 
resonant  wavenumbers ,  so  that  we  shall  under -estimate  the  final  peak 
height.  This  is  shown  in  Figure  It,  where  the  peak  is  shown  growing  to 
about  56  nm.  Brueck  and  Ehrlich  measure  their  peaks  at  about  100  nm.  The 
discrepancy  here  gives  us  an  indication  of  the  accuracy  of  our 
approximations.  If  a  method  can  be  determined  for  calculating  the 
"baseline"  independently  of  our  present  model,  which  is  obviously  not 
valid  for  off-resonant  wavenumbers,  then  we  should  be  able  to  more 
accurately  predict  the  maximum  peak  height. 

Our  discussion  so  far  has  revolved  around  the  third  of  the  three 
questions  raised  at  the  beginning  of  this  article.  We  would  now  like  to 
briefly  discuss  the  other  two.  We  have  considered  several  possibilities 
for  the  mechanism  by  which  the  plasmon  increases  the  deposition  rate.  The 
first,  which  at  one  point  we  considered  most  likely,  is  that  the  plasmon 
field  induces  a  dipole  in  the  adatom/molecule,  which  in  turn  attracts  it 
to  the  surface.  Unfortunately,  the  growth  rate  resulting  from  this 
mechanism  is  about  15  orders  of  magnitude  too  slow,  and  hence  it  is 
discarded. 

A  second  possibility  is  that  a  charge  transfer  occurs  between  the 
surface  and  the  cadmium  atom,  and  hence  the  plasmon  interacts  with  an  ion 
and  not  a  dipole.  This  appears  to  yield  a  growth  rate  on  the  right  order 
of  magnitude,  but  it  remains  unclear  how  the  charge  transfer  takes  place. 
It  is  conceivable  that  the  plasmon  field  enhances  the  charge  transfer 
process  also,  in  which  case  a  rate  proportional  to  the  square  of  the 
intensity  is  not  impossible. 

Yet  another  possibility  is  that  the  existence  of  the  plasmon 
radically  alters  the  shape  of  physisorpt ion/chemisorption  potential  well. 
This  raises  the  possibility  that  the  growth  is  not  proportional  to  the 
electric  field  intensity,  but  rather  to  the  electron  distribution  along 
the  surface.  The  very  complexity  of  this  problem  forces  us  to  exclude  it 
from  immediate  consideration. 

Finally,  there  is  the  dynamics  of  the  dissociation  process.  We  have 
reason  t^believe  that  this  occurs  in  the  gas  phase,  well  above  the 
surface.  If  this  is  indeed  true,  then  the  problem  becomes  considerably 
simpler  since  the  rate  would  depend  only  on  the  partial  pressure  of 
cadmium  in  the  gas  phase.  We  could  then  separate  the  problem  into  two 
parts:  the  gas-phase  problem,  i.e.,  the  production  of  cadmium  atoms,  and 
the  deposition  problem.  Depending  on  the  pressure  of  the  gas,  one  or  the 
other  will  predominate  in  the  total  rate  equation. 

This  concludes  our  discussion  of  recent  progress  made  in  the 
explanation  of  the  periodic  deposition  phenomena.  It  only  remains  to  make 
a  few  brief  comments  concerning  the  possible  applications  of  such  a  model. 
An  obvious  question  of  interest  would  be  how  one  could  maximize  either  the 
grating  height  or  growth  rate.  The  development  of  a  simple  model  will 
make  it  possible  to  predict  ways  to  do  this.  Similarly,  it  may  be 


possible  to  change  the  deposition  pattern  by  changing  the  parameters  of 
the  experiment.  The  calculation  of  this  phenomenon  involves  more  than 
what  we  have  discussed  here. 
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